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Abstract 

Bound (/(/-systems are considered in the framework of three different versions of the 
3-dimensional reduction of the Bethe-Salpeter equation, all having the correct one-body 
limit when one of the constituent quark masses tends to infinity, and in the framework of 
the Salpeter equation. The spin structure of confining qq interaction potential is taken in 
the form x'ji'y® + (1 — with < x < 1. The problem of existence (nonexistence) of 

stable solutions of 3-dimensional relativistic equations for bound qq systems is studied for 
different values of x from this interval. Some other aspects of this problem are discussed. 
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As it is well-known, at present the spin (Lorentz) structure of the qq interaction is not 
established theoretically in QCD, a fundamental theory of strong interactions. Consequently, 
it is interesting to consider different possible choices for the spin structure, as it was done 
in Refs. UJ-Q where the bound qq systems in the framework of the Salpeter equation were 
investigated. Further, it is well known that the Salpeter equation is the simplest version of the 
3-dimensional (3D) reduction of the Bethe-Salpeter (BS) equation, the latter being believed to 
provide a natural basis to study bound qq systems in the framework of the Constituent Quark 
Model. Namely, the Salpeter equation is obtained from the BS equation in a straightforward 
way, when the kernel of the latter is assumed to have the instantaneous (static) form. However, 
in the instantaneous approximation there exist other possible versions |5|, U of the 3D reduction 
of the BS equation which, unlike the Salpeter equation, have a correct one-body limit when 
the mass of one of the constituents tends to infinity. These versions, which can be derived 
by choosing appropriate effective 3D-Green's function for two noninteracting fermions, will 
hereafter be referred to as the MW and CJ versions, respectively. Moreover, a new version of 
the effective propagator for two free scalar particles, guaranteeing the existence of the correct 
one-body limit for 3D-equations, was suggested in Ref. @. The effective 3D-Green's function 
for two noninteracting fermions can be constructed along the lines similar to Ref. 0, in a 
standard manner (see below). 

Taking into account the fact that the relativistic effects are important for qq systems con- 
taining two light quarks, as well as for heavy-light systems, it seems interesting to carry out the 
investigation of this sort of systems in the framework of the above mentioned 3D relativistic 
equations and study the dependence of the properties of the bound qq systems on the spin 
(Lorentz) structure of the confining part of gg-interaction. In the present report, we deal with 
this problem as concerned to the qq mass spectrum. 

The relativistic 3D equations for the wave function of bound qq systems, corresponding 
to the instantaneous (static) BS kernel (K(P;p,p') — > K s t(p,p )), for all versions considered 
below can be written in a common form (in the cm frame) 



where M is the bound system mass. The equal-time wave function $m(p) is related to the BS 
amplitude $p(j>) by 



$m(p) = G 0eff (M,p) / = V(p,p')]® M (p') 




(1) 




(2) 
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and the effective 3D Green's function of the two noninteracting quark system Goe// is defined 
as 

G 0eff {M,p) = [ ^ [G Qeff (M,p) = goeff(M,p)(^ + m x )(fa + m 2 )] (3) 

Here, G 0e ff is the effective two fermion free propagator and goeff is the effective propagator of 
two scalar particles. The operator (?oe// can be given in the form 

Goe ff (M,p) = E E ^Tm^ A - 1Q2) ® ^ P = W ^ 

ai =±a 2 =± a \ 1V1 iP) 

where the projection operators A^ 1 " 2 '' are defined as 

Air 2) (pi,P 2 ) = Ai ai) (Pi) ® A? 2) (P 2 ), Al ai) (p1) = { ~ + Hm 

hiipi) = aS + rriiji, Ui = (m- + p 2 ) 1/2 (5) 
and the functions D^ aia2 \M, p) and d(M,p) are given by 

D^(M,p)= (pi F V rf ( M ^) = 1 

c^i + cj 2 — (tti-C/i + a 2 -c/2j 

2 2 

E 1 + E 2 = M, E l -E 2 = l — 2 = b (MW version) (6) 
D^\M,p) = (E 1 + aiwi)(£; 2 + a 2 u 2 ) 
d(M,p) = 2(uji + lu 2 ), a = Ef-tof = [M 2 + b 2 - 2(cof + u$)]/4 (CJ version) (7) 

(see Ref. || where the mass spectrum of the bound qq systems was investigated in the frame- 
work of the MW and CJ versions of the relativistic 3D equations in the configuration space). 

Note that expression (7) was derived from Eq. (3) by using that for g 0e ff(M,p) determined 
from the dispersion relation. The same relation is satisfied by the expression of go e ff{M,p) 
suggested in Ref. R] (see formula (10) therein). According to the prescription given in Ref. 
|7|], particles in the intermediate states are allowed to go off shell proportionally to their mass, 
so that when one of the particles becomes infinitely massive, it automatically is kept on mass 
shell and the corresponding equation can be reduced to the one-body equation. Using this 
expression for g Qe ff(M,p) from Eq. (3) we obtain the expression for Goe// given by Eq. (4). 
(Note that formula (11) in Ref. |7j is not correct as it does not follow from formula (10)). So 
we obtain 



£)(«i«a)(M,p) = {Ei + a 1 co 1 )(E 2 + a 2 uo 2 ) - 

2 



^' + (Ei + aiui) - (E 2 + a 2 tu 2 ) 



d(M,p) = 2RB, R=(b 2 - 4y 2 a) 1/2 , B 
b = M + b y, y 



R-b 



R-b 



+ a, 



(8) 



mi + m 2 

This version will below be referred to as the MNK version. 

Using the properties of the projection operators A^ 1 " 2 '* and formulae (4-8), the following 
system of equations can be derived from Eq. (1) 



[ M - ( ax u x + a^ 2 )]$ir 2) (p) = 
A^\M,p)A^\p,-p) f -^lh° 2 V(P,p') E E *m A \p 



(27T) 



(9) 



where §M ia2 \p) = A^ ia2 \p, —p) &m(p) an d the functions A^ aia2 ^ are given by 



(0102) 



A (±±) = ±1, 



M 



loi + u; 2 



^(aiaa) 



A ( aitt2 ) = M + (a^i + a 2 a; 2 ) 
2(a;i + u 2 ) 



a[M + (ai^i + a 2 uo 2 )}- 



(MW version) (10) 
(CJ version) (11) 



-[M - (<X\UJ\ + a 2 uj 2 )] 



R-b 



R-b 



+ {Ei + aiui) - (E 2 + a 2 to> 2 ) 



(MNK version) (12) 



As to the Salpeter equation, it can be obtained from the MW version by putting A^ ±T ^ = 
and $g T) = 0. 

Further, we write the unknown function l^ 1 " 2 '' in Eq. (9) in the form analogous to that 
used in Ref. 0, where the bound qq systems were studied in the framework of the Salpeter 
equation 

1 \ / 1 



where 



r( Q i«2), 



X ( M ia2 \p) (13) 



N- 



oavip/ (ooi + «imi) / V -a 2 a 2 p/ (oo 2 + a 2 m 2 ) 
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\P) 



ij ^Nr , {p)Nr\p) 



V 2a;i 7 V 2lu 2 
Then, if the qq interaction potential operator V(p,p') is taken in the form 



(14) 



V(p,p') = ^Vogip-P) + [Xlh°2 + (1 - x)I X h]V c {p-p'), (0 < X < 1), (15) 

the following system of equations for the Pauli 2 <g> 2 wave functions Xm can De derived 

[M-(a 1 u 1 + a 2 u 2 )] X ^ a2 \p) = 



A^\M;p)Y: E f^^ M) ^f\^)x^V) (16) 



Q l =± a 2 =± 



where the effective qq interaction operator V^} 012 ' 011 ' (p, p , <j\, cr 2 ) is expressed via the poten- 
tials V og and V c and some functions, taking account of relativistic kinematics. On using the 
partial-wave expansion 

X ( M a2 \p)= E <n\LSJMj> R^?\p), (n = p/p) (17) 

LSJMj 

From Eq. (16) the system of integral equations for the radial wave functions R^sj^ip) can t> e 
obtained. 

At the first stage our investigation is aimed at comparative analysis of the different versions 
of the 3D relativistic equations as concerns the existence of stable solutions for different values 
of the vector-scalar mixing parameter x in Eq. (15). For this reason, in Eq. (15) we neglect 
the one-gluon exchange potential and take the confining potential V c (r) in the oscillator form 
used in Ref. [10], which is a simplified though justified version (for the light and light-heavy 



sectors) of a more general form used in Ref. @ . Namely, we take 

V c (r) = A ~a s (rnl 2 )(^r*-V 

m x m 2 2 12tt / Q 2 \ x 

P>i2 = > m 12 =m 1 + m 2 , a s {Q ) = — — — ln-^- (18) 

mi2 33 — 2nf\ A 2 J 

In the momentum space, the system of integral equations for the radial functions Rlsj (p) 
with the above potential is reduced to the system of second-order differential equations. The 
solution of this equation is written in the form similar to that given in Refs. |], [K| 

oo 

r { lsT\p) = e c L Q s^R { :r\p) (19) 

?1=0 

where R n L{p) are the well-known oscillator wave functions. Then, the system of equations for 
R^isj (p) i s reduced to the system of linear algebraic equations for the coefficients C L a gj 2 \ 

M Cfe 1 = £ £ H%Z&JM) (20) 

aj,a 2 Li b n 

Here it is necessary to stress that the matrix H explicitly depends (except the Salpeter 
version) on the meson mass M we are looking for. Concequently, the system of equations (20) 
is nonlinear in M. 
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By truncating the sum in (19) at some fixed value N max , the eigenvalues M and the cor- 
responding coefficients C^sjn can be determined from the system of algebraic equations with 
the dimension 4:(N max + 1) (for the Salpeter case we have 2(N max + 1) equations), provided 
the procedure converges with the increase of N max . If the procedure does not converge, we 
interpret this as absence of stable solutions to the initial equations. As it has been mentioned 
above, the system of equations (20) is nonlinear in M for the MW, CJ and MNK versions of 
the 3D equations. For the solution we use the iteration procedure: at the first step, in the r.h.s. 
of Eq. (20) determining M, we substitute the solution of the Salpeter equation and then find 
M by the procedure explained above. At the second step, we substitute the obtained solution 
into the r.h.s. of Eq. (20) and the iterations are continued until the result converges. 

In the present paper, we calculate the masses for the following qq systems with nonequal 
mass quarks: ds (%, 3 Si, 1 P l , 3 P , 3 ^i, 3 ^2, 1 D 2 , 3 D h 3 D 3 ), cu and cs (%, 3 S U 1 P 1 , 3 P 2 ), 
for which the values of the meson masses are known. In the calculations the following values 
of the parameters of the qq -interaction potential (18) were used ||, [UJ: c^o = 710 MeV, 
V = 525 MeV, A = 120 MeV, m u = m d = 280 MeV, m s = 400 MeV, m c = 1470 MeV. 

On the basis of calculation of mass spectra of the above qq systems we have arrived at the 
following conclusions: 

The stable solutions of the MW, CJ and MNK versions of the 3D relativistic equations 
always exist for x = 0. For x — 1 these solutions do not exist for the majority of states 
under consideration. For the Salpeter equation the situation is just opposite: for x — 1 stable 
solutions always exist whereas for x = the solutions do not exist for the majority of the states 
studied. This agrees with the results obtained earlier in Refs. for the qq systems of equal 

mass quarks from the light quark sector (u,d,s). Moreover, for the CJ and MNK versions stable 
solutions always exist for (0 < x < 0.5). As to the MW version, this sort of solutions exist 
only for qq systems with one heavy quark, whereas in order to provide the existence of stable 
solutions in the same interval of x for light qq systems it is necessary to accept a much smaller 
value for the confining potential strength parameter ujq (e.g. 450 MeV instead of 710 MeV). 
As to the interval (0.5 <x< 1), the existence (nonexistence) of stable solutions in the MW, 
CJ and MNK versions depends on the quark sector (light or heavy), quantum numbers of the 
qq system and on the value of the parameter ujq. For the case of the Salpeter equation, the 
situation again is opposite - stable solutions always exist in the interval (0.5 < x < 1) for qq 
systems with both quarks from light-quark sector @-{|], in the whole interval (0 <x< 1) for qq 



5 



systems with both quarks from heavy-quark sector (cc) || , or from heavy-light sector (present 
result). The existence of stable solutions of the relativistic equations under consideration is 
mainly related to the presence of the mixed (H — and — h) energy components of the wave 
functions in the equation for the qq bound state. 

To illustrate these conclusions, in Tables 1. a,b,c we give the results of numerical solution 
of the system of equations (20) for ds, cu and cs bound systems for the states 1 S , , 3 Si and 1 P 1 . 

Note that in order to obtain stable solutions of the 3D relativistic equations, it is sufficient 
to take N max = 4 — 7 in the series (19). This property is common for all 3D versions and meson 
states under consideration. 

A more detailed analysis (and the comparison with experiment) of the results in the fra- 
mework of the above considered versions of 3D equations with the confining potential (18) 
(including the regularization problem of the wave function normalization condition in the CJ 
and MNK versions), as well as the description of decay properties of pseudoscalar and vector 
mesons (P — > fiu, V — > e + e~), will be published separately. 
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Russian Foundation for Basic Research under contract 96-02-17435-a. 
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Table l.a The mass spectrum (in GeV) for (ds) system 



oiaies 


vei sions 




x — u.u 


X — U. 1 


X — U.O 


v n ec 

X — U.O 


v— n 7 

X — U. 1 


v— n q 
x — u.y 


v— 1 n 

X — l.U 




l\/T"\Ai r + 

MW 


++ 
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o n 
all 
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1 n*78 
l.U / o 
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l.loU 
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1 C 


1V11N IV 


i i 

1 — h 


n oqs 
U.yoo 


u.yoo 


1.U10 


1.U0U 


1 nofi 
i.uys 


1 1 Q1 

1.101 


1 1 AP. 
1 .140 
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an 


n QQD 

u.yyu 


1 nnEc 
1 .uuo 


l.UOO 


1 DfiQ 

i.uoy 


* 


* 


* 
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1 — h 


u.yor 
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++ 
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i i 
++ 


n nscc 
U.yoo 


1 m 7 
1.U1 / 


l.U/ 
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l.ZZU 


1 OQn 
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» 
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all 


l.Uoo 
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l.Uoo 
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i.iyz 


* 


* 
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++ 
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l.Uoo 
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1 117 
1.11/ 
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1.UU4 
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l.UZU 
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1.U04 
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l.uyu 


1 1 Qn 
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* 
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MW 
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1 1 nt? 
l.lUo 


1 1 O/l 
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1 1 /in 
1.14U 


* 




t 

L/J 
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1 O EC8 
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oil 
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1 Q /I O 
l.o4z 
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1.381 


1.419 


* 


* 
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++ 


1.235 


1.274 


1.343 


1.403 


1.454 


1.499 


1.519 




» 


all 


* 


* 


1.326 


1.398 


1.453 


1.499 


1.517 



MW+ - uj q = 450 MeV 

* - absence of stable solution 
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Table l.b The mass spectrum (in GeV) for (uc) system 
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X — U. 1 


x — u.o 
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oil 
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z.lzy 


9 1 £7 
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* 
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z.UzU 
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9 1 /in 
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9 1 QQ 
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z.UOo 
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9 nm 
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* 
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z.UOO 


9 DR8 
Z.UOO 


9 i nn 
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z.l i U 




» 


oll 
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MW 
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++ 
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9 Qflfl 
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t 
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z.zoo 


9 901 


9 QQQ 
Z.OOO 


9 QQ9 

z.ooz 
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z.4o ( 




» 


oil 

all 


O Q9Q 


9 Q/in 
z.o4U 


z.oOO 


9 Qn/i 
z.oy4 


9 /I OK 
Z.4z0 


9 A 

z.4oy 


9 /i 7n 
z.4 ( y 


1 P 
n 


1V11N iv 


i i 

1 — h 


z.zOl 


9 97/1 
Z.Z ( 4 


o qio 
z.Olo 


z.oOo 


9 Qn/i 
z.oy4 


9 A 97 
z.4z ( 


9 A A Q 
Z.44o 




» 


all 


2.318 


2.328 


2.348 


2.371 


2.397 


2.429 


2.448 




Sal 


++ 


2.254 


2.278 


2.323 


2.363 


2.401 


2.435 


2.451 




» 


all 


2.250 


2.276 


2.322 


2.363 


2.401 


2.435 


2.451 



* - absence of stable solution 
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Table l.c The mass spectrum (in GeV) for (cs) system 



States 


Versions 


OL\OLi 


x=0.0 


x=0.1 


x=0.3 


x=0.5 


x=0.7 


x=0.9 


x=1.0 




M w 


1 1 

+ + 


O 1 7/1 
Z.Li'L 


O 1 QQ 

z.loo 


O 01 P 
Z.ZLO 


O O/l 
Z.Z4Z 


O Ofi7 
Z.ZO ( 


o om 
Z.Zyl 


QflO 

Z.oUZ 




» 


oil 

all 


o onn 
Z.ZOi) 


Z.ZlO 


O OQ 1 

z.zdl 


O O/l Q 

z.z4o 


o o7n 
z.z /U 


* 


* 




n t 


i i 
++ 

o n 
all 


Z.lol 

o oi n 
Z.zlU 


O 1 Ofi 

z.iyo 

o 01 n 

z.ziy 


O 00£ 

Z.ZZO 

O OQ£ 

z.zoO 


O O^Q 

Z.Zdo 

O 0£7 

z.zO / 


O 070 

z.z t\) 

O OQ 1 

z.zol 


Z.oUO 


O Q 1 7 
Z.ol ( 

* 


1 c 


1V11N rv 


i i 

1 — V 


1 71 
Z. 1 rf 1 


o 1 qc; 
z.loO 


O OI 

z.zlz 


O OQQ 

z.zoo 


O OftO 

z.zoz 


o oqc; 
Z.ZoO 


O OOfi 

z.zyo 






an 


o ono 
z.zuz 


ODQ 

z.zuy 


z.zzo 


0/10. 
Z.Z40 


OPK 

z.zoo 


* 


* 




oai 


i i 

1 — h 


17/1 
Z. 1 1 4 


1 QQ 
Z. loo 


01 P 
Z.ZLO 


0/10 
Z.Z4Z 


OfiQ 
Z.ZOo 


0Q1 

z.zy i 


o qoo 

Z.OUZ 




» 


an 


1 70 
Z. 1 1 z 


1 Q7 
Z. lo ( 


01 P 
Z.ZLO 


0/10 
Z.Z4Z 


OP7 
Z.ZO I 


0Q1 

z.zy i 


Q.D1 
Z.OUl 




MW 


1 1 

++ 


O 1 7(? 
Z.l (0 


o i m 
z.lyl 


O 01 Q 

z.zlo 


o O/l c: 
Z.Z40 


o o7n 
z.z to 


O Ofl/I 

Z.zy4 


o one 
Z.oUO 






oil 

all 


o o 1 n 
Z.Z1U 


O 1 Q 

z.zlo 


O OQO 

z.zoz 


Z.Zdo 


O 07Q 

z.z / o 










i i 
++ 


O 1 QQ 

z.loo 


O 1 QQ 

z.lyo 


O 007 
Z.ZZ I 


z.zoo 


O OQO 

z.zoz 


O Q07 


O Q 1 O 

z.oiy 






oil 

all 


O 1 o 
Z.Z1Z 


o oon 
Z.zzU 


O OQO 

z.zoy 


Z.ZOi) 


O OQ/1 

z.zo4 


* 


* 


3 C 


M1N Iv 


i i 
++ 


O 1 7Q 

z.l 


O 1 Q7 

z.lo ( 


O OI P. 
Z.ZLO 


o O/l n 
Z.Z4U 


O OPK 

z.zoo 


O OQQ 
Z.Zoo 


o ooo 
z.zyy 






all 

an 


o on/i 

Z.ZU4 


011 
Z.Z11 


007 
Z.ZZ 1 


o o/i 

Z.Z40 


OfiQ 
Z.ZOo 


* 


* 




oai 


i i 

1 — h 


1 7fi 
Z. 1 I 


1 Q1 

z. iy i 


01 Q 
Z.Zlo 


0/1 ^ 
Z.Z40 


070 
Z.z 10 


OQ/I 

z.zy4 


o qoc; 

Z.oUO 




» 


oll 

ail 


O 1 7/1 
Z. 1 ( 4 


O 1 QO 

z.loy 


O 01 7 
Z.Z1 I 


O O/l A 
Z.Z l ± l ± 


O 070 

Z.Z < u 


O OQQ 

z.zyo 


Z.oUO 




MW 


i i 
++ 


O /I QQ 

z.4oo 


O A £0 


O /1D7 

ZA\) ( 


z.Ooo 


O KPT 

Z.oo < 


o Knn 
z.Oyy 


O PA A 

z.ol4 




» 


oil 

all 


O A 07 


o r^np 
Z.OUO 


z.Ozo 


Z.04 / 


O c;7Q 

Z.O (o 


o c;oq 
z.Oyo 


o pon 
Z.OZO 




n t 


i i 
++ 


o a c;o 
Z.40Z 


O A 7/1 
Z.4 1 4 


O K 1 P 

z.010 


z.OOO 


o sno 
Z.OyZ 


O POP 

Z.OZO 


O PA Q 

Z.04o 




» 


oil 

all 


o /i no 
z.4yy 


o n i 
Z.011 


z.ooy 


O KPA 
Z.004 


Z.OyO 


O fi07 
Z.OZ ( 


O PA K 

Z.040 


1 p 
n 


1V11N IV 


i i 

1 — V 


O A QO 

z.4oz 


o a c;o 

Z.40Z 


O A QO 

z.4oy 


O r:o/l 
Z.0Z4 


O KK7 
Z.OO 1 


o c;q7 
Z.Oo ( 


o pr\o 
Z.OUZ 




» 


all 


2.481 


2.491 


2.511 


2.534 


2.561 


2.587 


2.605 




Sal 


++ 


2.438 


2.459 


2.497 


2.533 


2.567 


2.599 


2.614 




» 


all 


2.434 


2.459 


2.496 


2.533 


2.567 


2.599 


2.614 



* - absence of stable solution 
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